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Abstract 

We investigate formation of an apparent horizon (AH) in high-energy particle collisions in four- 
and higher-dimensional general relativity, motivated by TeV-scale gravity scenarios. The goal is 
to estimate the prefactor in the geometric cross section formula for the black hole production. We 
numerically construct AHs on the future light cone of the collision plane. Since this slice lies to 
the future of the slice used previously by Eardley and Giddings (gr-qc/0201034) and by one of 
us and Nambu (gr-qc/0209003), we are able to improve the prefactor estimates. The black hole 
production cross section increases by 40-70% in the higher-dimensional cases, indicating larger 
black hole production rates in future-planned accelerators than previously estimated. We also 
determine the mass and the angular momentum of the final black hole state, as allowed by the 
area theorem. 
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I. INTRODUCTION 


Several scenarios in which the fundamental Planck energy could be O(TeV) have been 

tl Q 

proposed. In these scenarios, our space is a 3-brane in large |1|] or warped extra dimen- 
sion(s), and gauge particles and interactions are conhned on it. If this is the case, a black 
hole smaller than the extra-dimension size is well described as a U-dimensional black hole 
centered on the brane (where D is the total number of the large dimensions ^), and its 
gravitational radius is far larger than that of a usual black hole with the same mass. This 
implies that such black holes could be produced using future-planned accelerators, because 
the gravitational interaction becomes dominant in particle collisions above the TeV scale, 
and the black hole production cross section 


(Jbh ~ vr [rh{‘2fi)Y 


( 1 ) 


becomes sufficiently large. Here, /i is the energy of each incoming particle in the center- 
of-mass frame of the collision and r/j(2/r) is the gravitational radius of a U-dimensional 
Schwarzschild black hole of mass 2/i, given by 

167rGi5(2/i) 


r/j(2/r) = 


( 2 ) 


{D — 2)U£)_2 

where Gd is the D-dimensional gravitational constant, and Ud _2 is the {D — 2)-area of a 
unit sphere. 

The phenomenology of black hole production in accelerators was hrst discussed in P| Tor 
reviews, see P; for a related issue of black hole production in cosmic rays, see e.g. p). 
There are four stages in the time evolution of a produced black hole. The hrst one is horizon 
formation in the particle collision. Next the balding phase follows, in which classical emis¬ 
sion of gravitational waves occurs, and the produced black hole relaxes to a D-dimensional 
Kerr black hole, whose metric was found by Myers and Perry P|. The third stage is the 
evaporation phase, in which the black hole evaporates due to the Hawking radiation and 
superradiance^. The particles emitted in this process are observed as the signals in acceler¬ 
ators. As the black hole evaporates, its mass approaches the Planck mass. In this Planck 
phase, quantum gravity effects will become important. 


^ Thus^ = 4: + n, where n is the total number of large extra dimensions. 

^ See |3 and references therein for an interesting recent discussion of the role of superradiance. 
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left shock right shock 


FIG. 1: Schematic picture of the spacetime of colliding high-energy particles. 


The Planck phase may lead to a number of unexpected phenomena as predicted by 
string theory, non-commutative geometry, etc. On the other hand, the hrst three phases 
are well described by classical or semi-classical gravity. Quantitative predictions concerning 
these phases are important, since such predictions are necessary to test the validity of higher¬ 
dimensional general relativity. Furthermore, since quantum gravity effects will be observed 
by the difference from the semi-classical signals, precise prediction of the semi-classical 
signals is required. 

Related to the evaporation phase, several studies of the greybody factors of D-dimensional 


black holes are available 


^ (see also [l^ for related issues). On the other hand, it is also 
necessary to investigate the process of the black hole production and its relaxation, because 
the black hole production cross section is directly related to the black hole production rate 
in accelerators. Furthermore, the classical gravitational radiation determines the mass and 
angular momentum of the final state of the produced black hole, which provides the initial 
conditions for the evaporation process. Hence, the study of the high-energy two-particle 
system is an important problem. 

The high-energy two-particle system in four dimensions has been discussed to some ex¬ 
tent before the appearance of the brane world scenarios. The metric of one high-energy 
particle was obtained by Aichelburg and Sexl by boosting the Schwarzschild black hole 
to the speed of light with hxed energy /i. The gravitational wave emission in the axisym- 
metric system of two combined Aichelburg-Sexl particles was studied by D’Eath and 
D’Eath and Payne (summarized in Q]). A schematic picture of the spacetime with 
two Aichelburg-Sexl particles is shown in Fig. 1. Two particles collide at the speed of light. 
The gravitational field of each incoming particle is infinitely Lorentz-contracted and forms a 
shock wave. Except at the shock waves, the spacetime is flat before the collision (i.e., regions 
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I, II, and III). After the collision, the two shocks nonlinearly interact with each other, and 
the spacetime within the future lightcone of the collision (i.e., region IV) becomes highly 
curved. The ultimate goal would be to clarify the structure of region IV. If this is possible, 
the black hole production cross section and the gravitational waves emitted in the relaxation 
process could be determined. But this analysis is difficult because of the quite complicated 
structure of the gravitational field, and no one has succeeded in deriving the metric in region 
IV even numerically. 

Nonetheless, one can estimate the lower bound on the black hole production cross section 
only with the knowledge of regions I, II, and III. This can be done by finding an apparent 
horizon (AHb because the AH existence is a sufficient condition for the event horizon (EH) 


n 


formation 1^ (assuming the cosmic censorship 


ip H). 


As mentioned 


in 


see 


also |17|), Penrose (1974, unpublished) constructed an AH on the slice m = 0,^ < 0 and 
n = 0,M < 0 in the head-on collision case in four-dimensional spacetime. Because this 
AH has intrinsic geometry of combined two flat disks, it is often called the Penrose flat 


disks. Eardley and Giddings Q extended the AH solution of Penrose flat disks to positive 
impact parameters. They analytically derived the maximal impact parameter hmax for the 
AH formation in a grazing collision in the four-dimensional spacetime. Subsequently, one of 
us and Nambu extended this analysis to higher-dimensional spacetimes. The values of 
braax ln H-dimensioual cases were obtained numerically, and they can be well approximated 


by ~ 1.5 X 2 ^V,,(2/i). 

The AH method provides a lower bound on the true collision cross section. This lower 
bound depends on the slice used to determine the AH and becomes larger if a future slice 
is chosen. Indeed, the maximal impact parameter of the AH formation will be larger for 
such a slice (simply because it is possible that, for a given impact parameter, an AH has 
not yet formed on the old slice, while it forms by the time a later slice is reached). The 
lower bound would asymptotically approach the exact cross section as we move into the far 
future. Because of the monotonic growth, the further we move into the future, the smaller 
the difference between the true cross section and the estimate provided by the AH method 
would become. 


After the works of |l3,13|, one of us raised doubts in the validity of the setup of the 
high-energy two-particle system [l^, because of possible strong curvature effects in colliding 
shocks. However, this problem was later shown to be an artifact of the unphysical classical 


4 




point-particle limit: for a particle described by a small quantum wavepacket large curvatures 


do not arise |20| (see also j21|). Roughly, if a wave packet of Planck size 5z is taken, 
curvature remains sufficiently small, while corrections to the Aichelburg-Sexl geometry are 
0{5z/rh) 1. This argument justified the use of the Aichelburg-Sexl two-particle system 

to compute the black hole production cross section in elementary particle collisions. See also 


22j | , 23| for other characteristics of incoming particles that could affect black hole formation. 


In light of the above discussion, the purpose of this paper is as follows. In the previous 


analysis Ql, the AH was constructed on the union of the two incoming shocks u = 
0,u < 0 and u = 0 ,m < 0 (referred below as the old slice). However, it is clear from Fig. 1 
that this slice is not at all optimal in the sense that there exist other slices within regions I, 
H, HI, located to the future of the old slice. Motivated by this observation, we proceed with 
the AH analysis on the slice of the future light cone of the shock collision plane, given by 
the union of the outgoing shocks u = 0, u > 0 and u = 0, u > 0 (referred below as the new 
slice). This slice is optimal in the sense that it is the future-most slice that can be taken 
without the knowledge of region IV. By this analysis, we will improve the lower bound on 


the cross section of the black hole production. In addition, using the area theorem 


m, 


we 


will hnd restrictions on the mass M and the angular momentum J of the hnal state (i.e. 
the produced black hole after the balding phase). This part of the analysis is new compared 
to 18[, and provides indirect information about the spacetime structure of region IV of the 
Aichelburg-Sexl two-particle system. 

This paper is organized as follows. In the next section, we explain the system setup and 
derive the AH equation and the boundary conditions in the new slice. Then we present the 
analytic solution of the AH equation in the head-on collision case and explain the numerical 
method for the more physically important grazing collision case. In Sec. HI, we present our 
numerical results. We summarize the results for the maximal impact parameter and discuss 
the mass M and angular momentum J of the final state, as allowed by the area theorem. 
Sec. IV is devoted to the summary and discussion. 

Our new lower bounds on ubh, rnost precise to date, are summarized in Table H. 
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II. APPARENT HORIZONS IN THE HIGH-ENERGY PARTICLE SYSTEM 


A. System setup 


We begin by reviewing the Aichelburg-Sexl metric, describing the gravitational held of 
a high-energy particle. Following the analysis in 17], [l 8 [, we use the metric of a massless 
point particle of |ll[ |l^ Il2i 1^ that is obtained by boosting the Schwarzschild black hole 
to the speed of light with hxed energy /i = 7 M. The result is 


ds^ = —dudv + dr^ + r^dQ,]^_^ + ^{r)6{u)du‘^, 


(3) 


<h(r) = 


(4) 


—21ogr {D = 4), 

2/(T) {D > 5). 

Here we adopt tq = as the unit of length, which is close to r/i(2/i). 

The delta function in Eq. indicates that the coordinate system is discontinuous at u = 0 , 
and that a distributional Riemann curvature (i.e., a gravitational shock wave) is located 
there. The continuous coordinates are introduced by 


u = u, 


V = 


V — 2 \ogr9{u) + u6{u)/r‘^ 


(D = 4 ), 


r = r 1 — 


V + 29{u)/{D — 4)r 
u 


D-A 


+ u9{u)/G^-^ (A>>5), 


r.D-2 


9{u) 


(5) 

( 6 ) 

(7) 

( 8 ) 


where 0* are the coordinates on the {D — 3)-sphere. The metric becomes 




ds^ = —dudv + 


u 




dG + 


1 - 


u 


,D-2 


9{u] 


( 9 ) 


where 9{u) denotes the Heaviside step function. Note that u = is a coordinate singu¬ 
larity for M > 0, because the {D — 3)-sphere shrinks to zero size. Thus the h > 0 region is 
mapped onto the 0 < u < region by this coordinate transformation. 

By causality, we can construct the metric of a high-energy two-particle system in regions I, 
H, and HI by simply combining the metric of the left and the right particles, because there is 
no interaction before the collision. Figure 2 shows the schematic spacetime structure adding 
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FIG. 2: Schematic picture of the spacetime of colliding high-energy particles with (D — 3) dimen¬ 
sions suppressed. The schematic shape of AH on the new slice (u > 0, u = 0 and u > 0, u = 0) is 
shown by solid lines, while the AH on the old slice (tt < 0, u = 0 and u < 0, u = 0) is shown by 
dashed lines. Dotted lines indicate coordinate singularities. 


one dimension to Fig. 1. Our goal is to construct an AH on the new slice, i.e., on the union 
of the two null surfaces u = 0, n > 0 and m > 0, n = 0. By the left-right symmetry (we 
work in the center-of-mass frame), it is sufficient to consider the u > 0, n = 0 surface. We 
introduce a coordinate 0 such that the metric in region II is given by 


ds'^ = —dudv + 


u 


1 + {D - 


dr'^ -|- 7-2 


u 


„D-2 


[d(j)^ + sin^ 0(ir2|,_4) . (10) 


The radial coordinate r in region II is adapted to the left particle, which is thus located 
at r = 0. In these coordinates, the right particle will cross the transverse collision plane 
M = n = 0 at a point distance b from the origin, where b is the impact parameter. We will 
choose coordinate^ so that this point is r = 6, 0 = 0. This setup is identical to the one 


used 


in |l7| 


and 


B. AH equation and bonndary conditions 

The schematic shape of the AH on the new slice is also shown in Fig. 2. Because u = r^~‘^ 
is a coordinate singularity, we have two boundaries in this analysis: at u = v = 0 and 
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FIG. 3: Schematic picture of the boundaries and Cout- We should solve for h{r, (p) in the region 
surrounded by the two boundaries. 

Cout at M = = 0. We show the schematic shapes of Cin and Cout in Fig. 3. Between 

these bonndaries, the AH shape is specihed by an nnknown fnnction u = h{r,(p). The 
tangent vector of the nnll geodesic congrnence of the AH snrface can be fonnd in terms 
of this fnnction nsing metric m and is given by 

1-2 / , \ -2 


= - 
2 

k^ = 2 , 


l + (D-3)- 


h 


,D-2 


hi + r 


-2 


1 - 


h 


, 0-2 


k^ = 


l + (F>-3)-^ 


-2 


hy , 


k^ = r-2 ( 1 


h 


^D-2 


-2 




( 11 ) 

( 12 ) 

(13) 

(14) 


Imposing that this congrnence has zero expansion, we get the AH eqnation: 


h, 


(r^-2-h) \h^rr + {D-3)^ 


{D -2)h- {3/2)rh^r , (T> - 2)h - (l/2)rh 

It 7-» O . / 1 T 


r^-2 + {D- 3)h 

hi 


r ^ ^ + (D - 3)h] <j + (D - 4) cot (ph^^ + 


rD-2 _ ^ 

D-3 D-7 


^D-2 


+ {D-3)h r^-^-h 


= 0. (15) 

Now we consider the bonndary conditions. At Cin, the continnity of the AH reqnires 

h{r,(p) = 0. (16) 

Another condition comes from the continnity of the nnll tangent vector (np to a factor). 
This condition is eqnivalent to k‘^{x)k'’{x*) = k^{x)k'^{x*) or 


K + ^ L ^ L* = 16 , 


(17) 
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where x* denotes the point symmetric to x with respect to the center of Cin (i.e. the point 
r = &/2, (j) = 0). 

Now we turn to the boundary conditions at Cout- Because u = r^~‘^ is a coordinate 
singularity, Cout has to be located at some fixed unknown radius r = Vmax so that the AH is 
continuous. Hence we have 

h = r^ax: ( 18 ) 

on Cout- 

The last boundary condition follows by imposing the continuity of on Cout- For this 
purpose, we should translate into the (h, v, f, 0) coordinates. Using the fact that h 
behaves like 

^ ^max ^,r(0)(^ '^max') T '^max') T ' ' ' (19) 

in the neighborhood of Cout, we obtain 

r = 1, (20) 

k‘ = [hr -(D- 2)r°3 /f, (21) 

yU ^Pmax 

k' = -rJri’ + 2(D - 2)-'hr/F, (22) 

k* = [ft,, - (B - 2)r°3 "" (r - r„„)-‘ft.,*/B. (23) 

where 

F = (D - 2)-X + 'ity [hr -(D- 2)C-^]h%. (24) 

In order that k^ be continuous, should be constant for all 0: 

[A,, - (B - 2)riy] /F = B. (25) 

There are also two other conditions given by 

= BrfCax cos (j), (26) 

o 3-D 

k^ = -^sin0, (27) 

r 

where we have used the symmetry of this system. (These conditions are analogous to the 
conditions for smoothness of a non-axisymmetric surface written in cylindrical coordinates.) 
Here, B and B are related as 

H = 1(28) 
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by the null condition = 0. Using Eqs. (El, ESD, El and (El, we derive 



(29) 


This is the second boundary condition at Cout- Although we have not used Eqs. (OH|l and 
(El, we can easily check the consistency. The boundary condition (I29|l is also consistent 
with the AH equation eg. Substituting the series El into Eq. ina, the leading-order 
term is 



+ {D -4) cot (ph^rej) + 


By substituting Eq. (El into the left-hand side, we can conhrm that this equation is actually 
satished. 

In summary, there are two boundary conditions for each boundary: Eqs. El and El 


for Cin and Eqs. El and (01 for Cout- We should determine the shape of the boundary 


Cin and the values of Tmax and B, as well as the function h(r, 0 ), so as to be consistent with 
these four boundary conditions. 


The AH equation El is highly nonlinear and Ending analytic solutions for b 0 seems 



almost impossible even for D = A. This is in contrast with the old-slice case |17|, where 
the AH equation was given by the Laplace equation. But as a numerical problem, the new 
case is quite similar to the old case except that one boundary is added (the Cin boundary 
conditions are the same as in the old case). We can solve this problem by extending the 


numerical techniques developed in m as explained later. 

C. Head-on collision case 

In the head-on collision case, we can solve the AH equation analytically. In this case, the 
function h depends only on r and the boundaries Cin and Cout are given by r = and 
r = Tmax) respectively. The equation and the boundary conditions become 
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h(Tmin) 0 , 
h,r(j’min) 2 , 
h{rmax) = r, 
h,r{rraax) = 2(-D - 2)r, 


.D-2 
max f 


D-3 
max • 


(32) 

(33) 

(34) 

(35) 


In D = 4 case, the solution is given by 


h = 2r^ log r, 

^min 1 ) 


^ max — v C • 


(36) 

(37) 

(38) 


In D > 5 cases, the solutions are as follows: 

2 


h = 


^D-2 . D-4 


T max — 


(D-4) 

= 1 , 

D _ 2 \ 


1 )> 


■) 


(39) 

(40) 

(41) 


(42) 


The total AH area is easily calculated. Restoring the length units, we have 

- Xd^) 

This is exactly the same value as the area of two Penrose flat disks (i.e., the AH on the old 

n 

slice) given in [121 • This coincidence can be interpreted as follows. Because the null geodesic 
congruence of the Penrose flat disk does not have shear, the expansion rate does not change 
(i.e. stays equals to zero) while propagating into region H according to Raychaudhuri’s 
equation. Hence, the null geodesic congruence of the AH on the new slice n = 0, m > 0 is 
the same as that of the Penrose flat disk, and the areas of the two AHs coincide. 

However, for positive impact parameters (grazing collisions), the shear of the null geodesic 
congruence of the AH on the old slice will be non-zero. While propagating, the expansion rate 
becomes negative according to Raychaudhuri’s equation, and the null geodesic congruence 
of the AH on the new slice will not be the same as on the old slice. This suggests that. 


using the new slice, we should hnd larger Ah 
compared to those in the previous results of 


areas and larger maximal impact parameters 


3,0 
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D. Numerical method for grazing collision 


In the grazing collision case, a numerical calculation is required to determine the AH. To 
solve the AH equation, we introduce coordinates (r, 0) by 

r cos(^ = r COS0 + 6/2, (43) 

r sin (j) = r sin (j). (44) 


In these coordinates, the central point of Cin is given by r = 0. (Because of the left-right 
symmetry, Cin will be symmetric with respect to the lines 0 = 0 and 0 = vr/2.) We specify 
the boundaries Cin and Cout by 


r = ^m(0), 
r = gouM) 


(1/2)6 cos 0 -f - (1/4)6^ sin^ 0 


1/2 


(45) 

(46) 


respectively. We make a further transformation 

R = , , (47) 

9out{4^') 9in ( 0 ) 

and solve the AH equation using the {R, 0) coordinates. The advantage of these coordinates 
is that Cin and Cout are specihed by i? = 0 and R = 1, respectively, and the boundary 
conditions are easily imposed. 

The following algorithm was used in the numerical solution of our boundary value prob¬ 
lem. First, some test boundaries Cin and Cout are taken, and the solution of the AH equa¬ 
tion satisfying only the two boundary conditions (USD and (HHD is constructed (using a 
hnite-difference discretization of the partial differential equation, and a convergent iterative 
procedure). Next, the difference from the boundary condition fITTjl is calculated: 

Am( 0 ) = {h‘^r + L L* “ ^ 6 , (48) 


and Cin is modihed as follows: 


= 9in(<P) + drAinW- 


(49) 


The Cout is also modihed at this step, as follows. Recall that Cout is characterized by Tmax 
and B. We determine B using Eq. (Ei at 0 = 0 and calculate the difference from the 
boundary condition (P|) at 0 = tt: 

= h,, - 2{D - 2)r^V(l - B). (50) 


12 





TABLE I: The error estimated by the difference from the Cout boundary condition, i.e., Eq. 
evaluated at 6 = bmax- The error decreases by a factor of about 4 when doubling the resolution. 
This indicates the correctness of our numerical calculations. 


D 

4 

5 

6 

7 

8 

9 

10 

11 

(25 X 50) grids 

0.21% 

2.0% 

4.2% 

- 

- 

- 

- 

- 

(50 X 100) grids 

0.052% 

0.49% 

1.1% 

1.9% 

3.0% 

4.0% 

4.6% 

5.1% 

(100 X 200) grids 

- 

- 

0.29% 

0.50% 

0.75% 

1.1% 

1.4% 

1.8% 


Using this value, we modify r^ax as follows: 


next 
' max 


= Vr, 


^out^out' 


(51) 


If we choose and eout appropriately, we can make the boundaries converge by iterating 
these steps. We truncated the iteration when the absolute values of and Aout become 

less than 10“®. 

To evaluate the numerical error, the following method was used. In the numerical method 
explained above, we could use only two grid points to impose the boundary condition at 
Cout^ because the values to be determined are only r^ax and B. In principle, the boundary 
conditions (El should be satished at the remaining grid points because Eq. El is compatible 
with the AH equation in the neighborhood of Cout, as expressed by El. However, in 
practice, because of the hniteness of the number of grid points, small differences from the 
boundary condition (El will be present in the intermediate grid points (0 < 0 < tt). This 
suggests to estimate the characteristic error as follows: 


5 = — V 

N ^ 


h^{D-2)r^J ( 1 + 


,r 


l-B^ 


1 + + 2B cos I 


(52) 


where the sum is taken over all grid points on Cout, and N is the total number of these 
points. 

Table I summarizes the resolution of the grids used to discretize the {R, 0) coordinates in 
our computations, as well as the error 6 a.t b = bmax- We observed that the error estimated 
by S becomes larger as b increases, and takes the largest value at 6 = bmax', it also becomes 
larger for larger D. For hxed b and D, the 6 typically decreases by a factor of about 4 if 
a grid with double resolution (i.e. with 4 times as many points) is used. Such behavior of 
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&= 0.4 



£>= 0.6 



6 = 0.8 



6 = 0.841 



-1.5 

-1.5 -1 -0.5 0 0.5 1 1.5 


FIG. 4: The shape of Cin and Cout for b/rQ = 0.4,0.6,0.8,0.841 in the D 
figure indicate the location of incoming particles. The old-slice AHs in 
are shown by gray lines. We could not find the AH for 6/ro = 0.842. 


= 4 case. Two dots in each 



for b/ro = 0.4,0.6, 0.8 


the error strongly indicates the correctness of our numerical program, and the convergence 
to the continuum limit. Although 6 seems somewhat large for large D, it turned out that 
the error in bmax is much smaller, as follows by comparing the values of bmax obtained for 
different grid resolutions. We estimate the error in bmax at the level of about 0.2% for all 
D. Instead, the 6 reflects the error in the AH shape. To summarize, the error in Figs. EllHl 
MM is roughly of the magnitude given in Table I, while the error in Table II is about 0.2%. 


III. NUMERICAL RESULTS 

In this section, we present the numerical results for the AHs in the grazing collision. The 
section is divided into two parts. We hrst provide the results for the AH shape and the 
maximal impact parameter of the AH formation. Next, we introduce two quantities that 
indicate the amount of energy trapped by the AH and discuss the hnal state of the produced 
black hole, as allowed or prohibited by the area theorem. 


A. AH shape and the maximal impact parameter 


Figure lU shows the shapes of Cin and (Tw for various values of b in the D = A case. 
The old-slice AHs of Eardley and Giddings [l7| are also shown. As b increases, becomes 
oblate, and Vmax becomes smaller. For small b, and the corresponding boundary curve 
of the old-slice AH almost coincide, which again indicates the correctness of our numerical 
program. For larger b, Cin lies outside the old-slice AH curve. For even larger b, we have 
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/j=0.4 ij=0.9 ^=1.1 fc=1.145 






FIG. 5: The shapes of Cin and Cout for b = 0.4,0.9,1.1,1.145 in the D = 5 case. Two dots in each 
figure indicate the location of the incoming particles. The old-slice AHs of [l^ for b/vQ = 0.4, 0.9 
are shown by gray lines. We could not find the AH for b/ro = 1.146. 


6 - 0.4 6 - 1.0 6 - 1.3 6 - 1.333 






FIG. 6: The shape of Cin and Cout for b/r^ = 0.4,1.0,1.3,1.333 in the D = 6 case. Two dots in 
each hgure indicate the location of the incoming particles. The old-slice AHs of jl^ for b = 0.4,1.0 
are shown by gray lines. We could not find the AH for bjro = 1.334. 


a situation when there exists an AH in the new slice, while there is no AH in the old slice. 
The become, about 3% latget thau the ptevioue teeult of Q. 

Figure El shows the shapes of Cin and Cout in the D = 5 case. In this case, Tmax is almost 
constant for all b and the shape of Cin becomes oblate as b increases. It is quite interesting 
that Cin becomes non-convex around b = bmnx- The value of b^ax is about 18% larger than 
the previous result of Yoshino and Nambu [ig, which leads to 40% larger cross section of 
the AH formation, the present value being ctah — I-Stt [r/i(2/i)]^. 

Figure El shows the shapes of Cin and Cout in the D = 6 case. In this case, r^ax becomes 
larger as b increases. The shape of Cin at 6 = bmax is even more non-convex than that 
in the D = 5 case. The value of bmax is about 26% larger than the previous result of 
ji^ . This leads to 59% larger cross section of the AH formation, the present value being 
(Tah ^ 2.171 [rh{2fi)f. 
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D=6 



FIG. 7; Three-dimensional plots of the AH shape function at 6 = bmax for D = 4, 5, and 6. The 
AH becomes taller for larger D as a consequence of the growing power exponent in the boundary 
condition m- 



b/ro 

FIG. 8: The relation between b/ro and where rmin = The maximal impact 

parameter occurs at drmin/db = —oo. 

In Fig. 0 we plot the AH shape function h(r, 0) at 6 = bmax for D = 4,5, and 6. The AH 
becomes taller for larger D as a consequence of the growing power exponent in the boundary 
condition (HED- 

For D > 7, the shapes of Cm and Cout and the horizon shape behave qualitatively similarly 
to the D = 6 case, and we do not present them here in detail. 
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TABLE II: The value of ftmax/^’O; the ratio of the increase in the maximal impact parameter, and 


the value of the AH formation cross section ctah (which provides a rigorous lower bound for (Tbh)- 


D 

4 

5 

6 

7 

8 

9 

10 

11 

bmax ! “^0 

0.841 

1.145 

1.333 

1.441 

1.515 

1.570 

1.613 

1.648 

bmax ! bmax 1 

5% 

18% 

26% 

29% 

29% 

30% 

30% 

30% 

o-ah/tt K(2p)]^ 

0.71 

1.54 

2.15 

2.52 

2.77 

2.95 

3.09 

3.20 


^max/^h^ /^) 

1.75 
1.5 
1.25 
1 

0.75 
0.5 
0.25 


4 5 6 7 8 9 10 11 



FIG. 9; Summary of the b^nax/'f’h{‘^^^) values (o) and the previous values (x) for 

17 = 4,..., 11. 

Figure |H1 shows the minimum radius of Cin, Tmin = l‘^)) as a function of h for 

D = A ,..., 11. The maximal impact parameter occurs at drmin/db = —oo. From this hgure, 
we can read off the values of bmax- 

Table II summarizes the numerical results for bmax, which are the most important output 
of our analysis. For 5 < Z7 < 11, the values of b^ax increase by 18-30% compared to 
the previous values of bmax in U. Con.espo„di„g.y, tKe value, of the etos. seet.ou of 
the AH formation ctah/tt [rft(2p,)]^ increase by 40-70 %. This indicates that the black hole 
production rate in accelerators can be quite a bit larger compared to the previous estimates, 
this tendency being especially enhanced for larger D. In the D = A case, which may have only 
astrophysical applications, we hnd only a modest 5% improvement in bmax compared to |l7| . 
We compare the present values of bmax/^hi^^l^) with the previous values of bmax/^hi^^l^) in 
Figure El 
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B. Trapped energy and final state of the produced black hole 


found AHs (see 


0 


lip 0. 


Assuming the cosmic censorship an event horizon (EH) must be present outside the 


for a proof). Moreover, the area theorem states that the EH area 


never decreases. Hence one naturally expects that the AH mass dehned by 


Mah — 


{D — 2)flD-2 f Ad- 


IbvrG 


D 


iD-2 

^D-2 


{D-3)/{D-2) 


(63) 


provides the lower bound of the irreducible mass Mirr of the hnal Kerr black hole. We 
should mention that for this statement to be rigorously justihed, the area of an arbitrary 
surface outside of the AH in a given slice should be larger than that of the AH. In the old 


slice analyzed in HQ this property actually holds. On the other hand, on the new slice 
analyzed here we can End a counterexample. In the head-on collision case, for example, the 
union of two surfaces r = Tc > Tmax, 0 < u < u = 0 and r = rc,0 < v < u = 0 is 
closed, lies outside of the AH, and has zero area. Thus there is no rigorous reason why in 
the present analysis Mah should provide the lower bound on the final mass. 

Nevertheless, we can still hnd a rigorous lower bound on Mj^r, arguing as follows. The 
intersection of the AH and u = u = 0 is given by Cin- Let us denote the intersection of the 
EH and u = v = 0 hy E'er- This curve must lie outside Further, one can show that the 
area Aeh of the intersection of the EH with the old slice is equal to twice the area of the 
region surrounded by Ceh calculated with the {D — 2)-dimensional flat metric. It follows 
that Aeh is bounded below by Aib, the latter quantity being dehned as twice the area of the 
region surrounded by Cin, calculated with the same hat metric. Hence the rigorous lower 
bound on M^rr is given by 

(D-3)/{D-2) 


Mib = 


{D — 2)f2E)-2 


IbvrG 


D 


A, 


lb 


^D-2 


(54) 


This Mib as a function of b is shown in Fig. [THl The non-rigorous bound Mah is also 
included, because under some additional assumptions it may still provide the energy trapped 


by the produced black hole. For example, the Hawking quasi-local mass 


26| calculated on 


the AH coincides with Mah- The old-slice value of the AH mass Mah as found in |17|, ll8| is 
also shown. We see that Mib and Mah fake close values, Mah being slightly larger. Although 
Mib and Mah are close to Mah for small b, they becomes significantly larger around b = bmax, 
especially in the higher-dimensional cases. 
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FIG. 10; The rigorous lower bound Mu, of the final irreducible mass (black lines) and the indicator 
of the trapped energy Mah (light gray lines) for D = 4,11. The previous value of the AH mass 
Mar of 12, 3 are also shown by dark gray lines. 


Now we consider the mass M and the angular momentum J of the hnal Kerr black hole 
which are allowed by the area theorem: 


Mirr > Mil,. ( 55 ) 

Here = Mirr{M, J) is the irreducible mass of the Kerr black hole, which is defined, just 
as in four dimensions, as the mass of a Schwarzschild black hole having the same horizon 
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area. It is thus related to the Kerr black hole horizon area A^err by the formula: 


Axerr — ^{Atirr)- 

The left-hand side of this equation can be easily computed using the explicit D-dimensional 
Kerr black hole metric [3], which gives the relation 


r^-\M,rr) = r^-\M)n{M, J). (56) 

Here rk{M, J) is the Kerr black hole horizon radius which satisfies the following equation 


rl{M, J) + 


{D - 2)J 
TM 




5-D/ 


(57) 


The total energy and the angular momentum of the system before the collision are 2/i 
and &/i, respectively. Denoting 


e = M/2p, 

(68) 

C = 

(59) 


the hnal state of the produced black hole will be specified by a point in (,^, (C)-plane with 

0 < ^ < 1 and 0 < C < 1. 

For D = 4, 5 there exists an upper limit on the black hole angular momentum for a fixed 


J < MM) 


(l/2)Mr,(M) (D = 4), 
(2/3)Mr,(M) {D = 5), 


(60) 


Region (i) consisting of points satisfying the opposite condition J > MM) should be a 
priori excluded from the (.^, (C)-diagram. 

According to the above discussion, region (ii) corresponding to black holes violating the 
area theorem (ESD can also be excluded. The remaining points constitute the allowed region 
(hi). 

Figures HTHTH show regions (i), (ii), (hi) for D = 4, 5, 6, and 9 for some selected values of 
h. We see that the condition > Mu, gives a stronger restriction on the hnal state (^, C) 
than the simple condition M > Mu,. This difference becomes quite noticeable especially for 
b ~ bmax in the D = 4 and 5 cases. This result indicates that M should be quite a bit larger 
than Mih at 6 ~ b^ax, because almost 100% angular momentum should be radiated away if 
M ~ M;fe, which would be quite unnatural. 
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FIG. 11: The regions (i), (ii) and (iii) in the (^,C)-plaiie for b = 0.4,0.6,0.8 in the D = 4 case. 
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FIG. 12: The regions (i), (ii) and (iii) in the (^, C)-plane for b = 0.5,1.0,1.145 in the D = 5 case. 

Unfortunately, we cannot find a non-trivial upper bound for the angular momentum J of 
the hnal Kerr black hole. (If the boundary of region (iii) intersected the ^ = 1 line at ^ < 1, 
we would be able to hnd such a bound.) On the other hand, it is interesting to note that our 
results are quite consistent with previous numerical simulations of gravitational collapse of 
rapidly rotating bodies in four dimensions (see 27| and references therein). In these works, 
the authors found a necessary condition for black hole formation expressed as 


Q — '-System/ system) ^ 1; 


( 61 ) 


where Mgystem and Jsystem are the total gravitational mass and angular momentum of the 
system. In our system, the value of g at 6 = hmax is 


q = 


(62) 


0.84 {D = 4), 

0.93 {D = 5), 

which is in agreement with It should be pointed out that for the hve-dimensional black 
ring solutions there is no upper bound on q, and thus we expect that criterion p[|) in 
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FIG. 13: The regions (ii) and (iii) in the (^, C)-plaiie for b = 0.5,1.0,1.3 in the D = 6 case. 
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FIG. 14: The regions (ii) and (iii) in the (^, C)-plane for b = 0.5,1.0,1.5 in the D = 9 case, 
the five-dimensional case holds only for formation of the AH with spherical topology. 


IV. SUMMARY AND DISCUSSION 


In this paper, we have analyzed the AH formation in the high-energy particle collision 
using a new slice u = 0 ,u > 0 and u = 0 ,m > 0 , which lies to the future of the slice 

nn 

u = 0,v < 0 and u = 0 ,m < 0 used in the previous studies of jl2l, I1SI|- Our main results 
are summarized in Table II. Compared to the previous results for bmax, we have obtained 
maximal impact parameters b^ax of the AH formation larger by 18-30% in the higher¬ 
dimensional cases. These results lead to 40-70% larger cross section of the AH formation, 
the present value being ctah — Stt [r/i(2p)]^ for large D. 

We have also estimated the mass M and angular momentum J of the hnal state of the 
produced black hole, as allowed by the area theorem Mirr > Mu,. This condition provides a 
stricter restriction on the hnal M and J than the simple condition M > Mu,, and becomes 
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especially effective for large b in the D = A and 5 cases, when onr results indicate that the 
hnal mass M should be signihcantly larger than Mu,. We also found that Eq. (El gives a 
necessary condition for the AH formation in the D = A and 5 cases, which is consistent with 
the various numerical simulations of the gravitational collapse. 

Our analysis provides the most precise data on the cross section of the black hole produc¬ 
tion in high-energy particle collisions to date. Using our new results, various phenomeno¬ 
logical discussions that used the results of |l^ (such as e.g. or relied on the more rough 
estimate m (e.g. Q and many others) can be improved. The present investigation is a 
necessary step towards the hnal understanding of the semi-classical signals that would be 
observed in future-planned accelerators. 

It should be stressed that the estimates on M and Mirr provided by our analysis give only 
rigorous upper bounds on the amount of emitted gravitational radiation. The real amount 
is likely to be smaller than suggested by these estimates, by a factor of a few. The work 
of D’Eath and Payne gives an idea about the size of this effect. In their analysis of 
axisymmetric collision of two Aichelburg-Sexl particles, they calculated the evolution of the 
gravitational held far away from the center using 7 “^ as a small parameter, and derived the 
news function near the symmetry axis to the second order. Assuming the azimuthal pattern 
of the gravitational radiation, they estimated the energy loss to be 16%, which should be 
compared to the rigorous upper bound of 29% provided by the AH method^. It is natural to 
expect that reduction of comparable size will occur in all dimensions. It should be mentioned, 
however, that a recent calculation based on an “instantaneous collision” approximation 
in linearized gravity predicts that gravitational wave emission becomes highly suppressed in 
higher dimensions (up to 0.001% in D = 10), which in our opinion is unlikely. Still another 
setup models the collision by a lightlike particle falling into a Schwarzschild black hole 
and gives estimates which are closer to our values ( 8 % in D = 10). We point out that all 
these works have problems such as ignoring the nonlinearity of the system, or the setup is 
too far from the realistic one. Analysis without approximations remains an important open 
problem. 

The ultimate goal of such analyses, which is left for the future, is to determine the 


^ It should be also noted that D’Eath and Payne’s estimate does not take into account additional gravita¬ 
tional radiation from the center of the system, which cannot be evaluated by this method. 
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spacetime structure after the collision, i.e., in region IV (u > 0,n > 0). This would clarify 
the precise maximal impact parameter of the black hole formation and the relation between 
the values of (M, J) of the hnal state and the impact parameter b. If this is completed, we 
will be able to obtain quite accurate semi-classical predictions by using the existing studies 
of the greybody factors. 
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